Controlling Viral Capsid Assembly with Templating 
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We develop coarse-grained models that describe the dynamic encapsidation of functionalized nanoparticles by 
viral capsid proteins. We find that some forms of cooperative interactions between protein subunits and nanopar- 
ticles can dramatically enhance rates and robustness of assembly, as compared to the spontaneous assembly of 
subunits into empty capsids. For large core-subunit interactions, subunits adsorb onto core surfaces en masse 
in a disordered manner, and then undergo a cooperative rearrangement into an ordered capsid structure. These 
assembly pathways are unlike any identified for empty capsid formation. Our models can be directly applied 
to recent experiments in which viral capsid proteins assemble around the functionalized inorganic nanoparti- 
cles [Sun et al., Proc. Natl. Acad. Sci (2007) 104, 1354]. In addition, we discuss broader implications for 
understanding the dynamic encapsidation of single-stranded genomic molecules during viral replication and for 
developing multicomponent nanostructured materials. 



I. INTRODUCTION 

The assembly of simple building blocks into larger, ordered 
structures is crucial for many biological processes and is en- 
abling novel nanostructured materials (e.g. |H IH SI [H |(| 
7]), which often draw inspiration and materials from biology. 
The spontaneous assembly of viral proteins into empty cap- 
sids (protein shells) has been the subject of elegant in vitro 
experiments (e.g. |§ i M O [IlJllJMJlslJlJJ^ and 
insightful theoretical works (e.g. HIMMlaHllll 

s s m m in, m Eh). The in vivo replication of many 

viruses, however, involves simultaneous assembly and encap- 
sidation of the viral genome ll32h . Likewise, many nanostruc- 
tured materials require precise spatial ordering of multiple, 
dissimilar components. In this article, we develop coarse- 
grained models for a particular example of multicomponent 
assembly - the assembly of viral capsid proteins around a 
rigid spherical template. Our models predict that a template 
enables increased assembly rates and efficient assembly over 
a wider range of parameters that control assembly driving 
forces, as compared to the spontaneous assembly of empty 
capsids. We find that template properties can control assem- 
bly pathways, and our models predict a novel mechanism that 
is unique to multicomponent assembly. 

Our models are motivated by recent experiments in which 
Brome mosaic virus (BMV) capsid proteins dynamically en- 
capsidate functionalized inorganic nanoparticle cores, cre- 
ating unique biological and synthetic composite structures 
called virus-like particles (VLPs) 11331 l34i l35l 13611 . By com- 
bining the unparalleled self-assembly and targeting capabil- 
ities of viruses with the functionalizability of n anop articles, 
VLPs show promise as imaging agent s If36l r3~7l l38l f39ll . di- 
agnostic and therapeutic vectors Il40l |4U 14211 . and as sub- 
units or templates for synthesis of advanced nano-materials 
iSHQSHEl- Our models offer a framework with which to 
interpret experimental results in order to design more efficient 
templated assembly of nanomaterials, and a means to use this 
as a model system with which to understand aspects of viral 
protein assembly around nucleic acid cores. 

Formation of the hollow shell geometry of a capsid poses 
a significant challenge that requires anisotropic, directional 
interactions between subunits. Thus, in addition to their 



biomedical and technological applications, studying viral cap- 
sids has revealed fundamental principles of assembly. Al- 
though specific assembly mechanisms are poorly understood 
for most viruses, a general mechanism has emerged for the 
spontaneous assembly of empty capsids ifToL [TH UH EH EH 
Q3, E3. EH H3. S3. ES S3. Eai. Assembly occurs through 
a sequential addition process in which individual subunits or 
larger intermediates l26l S3l bind to a growing capsid. As- 
sembly rates must be restrained to avoid two forms of kinetic 
traps (long-lived metastable states): (a) if new intermediates 
form too rapidly, the pool of free subunits becomes depleted 
before most capsids finish assembling lfT2L I2TL I26L I47I I48L I49L 
15011 . (b) malformed structures result when additional subunits 
bind more rapidly than strained bonds can anneal within a 
partial capsid ll26l EM Bll l52h . The formation of too many 
partial capsids can be suppressed by a slow nucleation step 
04811 . but avoidance of both sources of kinetic frustration re- 
quires relatively weak subunit-subunit binding free energies 
$M HI |H SI SI HI • Theoretical work suggests that weak 
binding free energies are a general requirement for success- 
ful assembly into an ordered low free energy product; binding 
free energies that are large compared to the thermal energy 
(kBT) prevent the system from 'locally' equilibrating between 
different metastable configurations during assembly Il53ll54ll . 

Although in vitro studies of empty capsids provide a foun- 
dation for understanding assembly, interaction of proteins 
with a central 'core' is crucial for the replication of many 
viruses in their native environments, where capsid proteins 
must encapsidate the viral genome during assembly ll32ll . 
There is no role for exogenous species in the sequential as- 
sembly mechanism discussed above, but in vitro capsid as- 
sembly experiments in the presence of RNA demonstrate dif- 
ferent kinetics than capsid proteins alone, and suggest the 
presence of protein-RNA intermediates 15511 . 

Prior theoretical and computational studies of multicom- 
ponent assembly have examined the equilibrium behavior of 
polyelectrolyte encapsidation 11561 l57l l58l 15911 and the equi- 
librium configurations of colloids confined to convex surfaces 
1I27I, l60ll . A qualitative kinetic model has been proposed to ex- 
plain the formation of icosahedral symmetry in encapsidated 
RNA [61] and Hu and Shklovskii 16211 considered a model in 
which capsid proteins nonspecifically bind to single-stranded 
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RNA and slide on it towards an assembling capsid at one end, 
which increases the rate of assembly. While increased binding 
rates are one possible feature of multi-component assembly, 
an interior core, such as a nucleic acid or nanoparticle, may 
also promote assembly by acting as a template that steers as- 
sembly towards certain morphologies and as a heterogeneous 
nucleation site that localizes capsid proteins in an environment 
favorable for assembly. These factors may generate assem- 
bly mechanisms that are entirely different from the sequen- 
tial mechanism considered in the formation of empty capsids. 
For instance, McPherson 1 63] proposed a qualitative model in 
which a large number of proteins non-specifically bind to a 
nucleic acid molecule to form a structure resembling a reverse 
micelle, and then reorient to form an ordered capsid. In this 
work we present a computational model for the encapsidation 
of an interior core with no pre-assumed pathways. For some 
sets of system parameters, our simulations predict assembly 
mechanisms consistent with McPherson's model. 



In addition to the technological applications discussed 
above, solid nanoparticles offer a simplified, controllable ex- 
perimental system with which to test models for the effect 
of heterogeneous nucleation and templating on assembly, and 
thus may lead to valuable insights about viral assembly around 
nucleic acid cores as well as elucidate the fundamental prin- 
ciples of multicomponent assembly. Experiments show that 
capsid assembly around cores competes with spontaneous as- 
sembly at subunit concentrations well above the threshold 
concentration for empty capsid assembly (critical subunit con- 
centration, CSC) 113311 . and that assembly occurs in the pres- 
ence of cores below the CSC 16411 . These results point to 
the ability of nanoparticles to act as heterogeneous nucleation 
agents. In addition, nanoparticles promote formation of cap- 
sid morphologies that are commensurate with nanoparticle 
sizes, which suggests that cores can direct the final assem- 
bly product through templating. The time dependence of the 
mass averaged amount of proteins on cores can be estimated 
by light scattering [64], but is not possible to characterize the 
extent to which these proteins have assembled without static 
procedures, such as crystallography or electron microscopy. 
Kinetic models that relate assembly pathways to dynamical 
observables such as light scattering are therefore necessary to 
understand assembly mechanisms. 



In this work, we present a computational model for assem- 
bly around solid cores, with which we analyze kinetics and 
assembly pathways as functions of the parameters that control 
the driving forces of assembly, including subunit concentra- 
tions, subunit-subunit binding energies, and surface adsorp- 
tion free energies. At low adsorption free energies and/or low 
subunit concentrations, assembly mechanisms resemble those 
seen for empty capsids, whereas assembly pathways at high 
adsorption free energies and/or subunit concentrations resem- 
ble the reverse micelle model. We demonstrate that the effect 
of cores on rates and assembly mechanisms can be understood 
through simple and general scaling arguments. 



n. MODEL 

We consider a dilute solution of capsid subunits with a re- 
duced concentration C s = pa 3 , with p the number density and 
a the subunit diameter, and rigid cores with a reduced con- 
centration Cq. Subunits can spontaneously assemble to form 
empty shells with well-defined structures of size N subunits. 
In addition, subunits interact favorably with cores and thus 
adsorb to, and assemble on, core surfaces. Complete assem- 
bly of adsorbed subunits results in core encapsidation. Our 
models are motivated by the experiments described above in 
which viral capsid proteins assemble on inorganic nanopar- 
ticles; thus, we begin by adapting a model previously used 
to simulate the spontaneous assembly of empty capsids j26Tl . 
This computational model is general, however, and could de- 
scribe, for example, colloidal subunits with directional inter- 
actions 1651 l66l 16711 . Likewise, the scaling arguments below 
are general enough to describe many forms of simultaneous 
assembly and cargo encapsidation for systems in which cargo 
degrees of freedom change slowly in comparison to assembly 
timescales. 



A. Modeling empty capsid formation 

We imagine integrating over microscopic degrees of free- 
dom as capsid proteins fluctuate about their native states, to 
arrive at a pairwise decomposable model in which subunits 
have spherically symmetric excluded volumes and direction- 
ally specific, short ranged attractions between complemen- 
tary interfaces. The lowest energy states in the model corre- 
spond to separate 'capsids', which consist of multiples of 60 
monomers in a shell with icosahedral symmetry. In this work 
we model experiments in which BMV capsid proteins assem- 
ble around 6 nm nanoparticles, for which only Tl capsid ge- 
ometries are observed 13311 . Because the basic assembly unit 
of BMV is believed to be a dimer and capsid proteins rapidly 
dimerize in solution 1 68], our model subunit represents a pro- 
tein dimer. Our energy minimum model capsid therefore is 
comprised of 60 monomers or A = 30 dimer subunits ar- 
ranged with icosahedral symmetry, as shown in Fig.Q] 

The locations of subunit interfaces are tracked by inter- 
nal bond vectors, h[ a \ that are fixed rigidly within a sub- 
unit frame of reference, with a 6 {1, 2, . . . , rib} an d n b is 
the number of interfaces on each subunit. In this work there 
are rib = 4 bond vectors that can be represented in Carte- 
sian coordinates as bj Q) /b = 0.5/c^x + 0.809fc£y + 0.309z 
with the bondlength b = 2~ 5/6 and k x = {1,-1,-1,1} and 
k y = {1, 1, —1, —1}, so that the angles between bond vectors 
have the values indicated in Fig.QJ and the minimum energy 
capsids have 30 subunits as shown in Fig.QJ). This model re- 
sults from merging pairs of monomeric subunits in the B5 cap- 
sid model considered in Ref.|26| (see Fig. 1 in that reference) 
and the resulting model capsid has the same connectivity as 
a model considered by Endres et al. 14811 (see Fig IB in that 
reference). 



The interaction between subunits i and j is 

Uij = u Kp QRi - R,-|) + «* (|rj a) - rf |) (1) 

a/3 

x s a p(e { ^\6 m ) Sa ^^\<i> m ) (2) 

where is the center of subunit i and r[ a ^ = Hi + b^ is 
the position of the interface represented by bond vector a on 
subunit i. The repulsive excluded volume interaction is 

u rep (i?) = 6(r m - R)[l + C(R/a)}, (3) 

where we have defined a 'Lennard- Jones' function C(x) = 
4(a; -12 — x~ 6 ) and 0(a;) is the step function with r rn = 2 1 / 6 <r 
the maximum range of the repulsion. 

The attractive interactions depend on the relative configura- 
tions and alignments of complementary interfaces. To reflect 
this, the sum in Eq. (f2]i is marked with a prime to note that 
it runs over all pairs of complementary interfaces, which in 
this work are all pairings of bond vectors with an even and 
odd label, e.g. (a, [3) = (1, 4) or (3, 2). These pairs are de- 
noted primary interactions. A favorable interaction between a 
pair of complementary interfaces has three requirements (see 
Fig- BO- First, the interfaces should closely approach each 
other, which is enforced by the distance potential 

u«t(r) = £b6(r c - f) [C(f/a) - £(r c /a)} (4) 

with r c = 2.5 the cutoff distance, £b the strength of the attrac- 
tive interaction, and f = r + r m is a shifted distance so that 
the minimum of the potential occurs at r = 0. The second 
requirement for a favorable interaction is that primary bond 
vectors,!^ and are aligned antiparallel, which is en- 
forced by the first factor on the second line of Eq. (0]i with 

cos(^)=-bf)-bf/6 2 - (5) 

The tolerance of the potential to angular fluctuations is con- 
trolled by the specificity function s a p: 

8 a p(i>, 1pm) = ^S(?/>m - VO[cOs(7n/V-0m) + 1] (6) 

where ?/; max is the maximum angle deviation, which controls 
the angular specificity of the attractive interactions. In this 
work, we vary the primary bond angle tolerance, 9 m , from 
0.25 to 2.5 radians. 

The third requirement for an attractive interaction is that 
two secondary bond vectors, which are not involved in the 
primary interaction, are coplanar. This requirement enforces 
angular specificity in the direction azimuthal to the primary 
bond vectors. A specific pair of secondary bond vectors (7, e) 
is associated with each primary pair (a, (3), and the second 
angular factor in Eq. (f2]i favors the alignment of the normals 
to two planes. The first plane is defined by the inter-subunit 
vector R.y = R; — R 3 and the first member of the secondary 

pair, b' 7 '' ; the second plane is defined by R^ and the second 
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FIG. 1 : (a) Geometry of subunits and attractive interactions for the 
computational model. Bond vectors are depicted as arrows; bond 
vector 3 on subunit A and bond vector 2 on subunit B have an attrac- 
tive interaction if they are nearly antiparallel and the secondary bond 
vectors, bond vector 2 on subunit A and bond vector 3 on subunit B, 
are nearly co-planar. The angles between indicated bond vectors are 
specified in degrees, (b) The low free energy capsid geometry. The 
center of each subunit roughly corresponds to a 2-fold axis of sym- 
metry (at dimer interfaces) in a BMV Tl capsid (see Ref. 71 and the 
VIPER database 1 72]). Subunit sizes are reduced to aid visibility. 

member of the same secondary pair, hj . Denoting these nor- 
mals by n] = b- 7 ^ x R^ and nj = x Ry , the dihedral 
angle <j> in Eq. (O is determined from 

cos(^f) = n! 7) -nf (7) 

with n the unit vector n/|n|. 

Subunit positions and orientations are propagated accord- 
ing to overdamped Brownian dynamics, with the unit of time 
to = a 2 /48D, where D is the subunit diffusion coefficient. 
All energies are measured in units of the thermal energy k%T. 
In this work, secondary pairs are specified as the inverse of 
the primary pairs: (7, e) = ((3, a), and the dihedral specificity 
parameter <fi m — it radians throughout. We vary the subunit 
concentration over the range Cs G [2 x 10~ 3 ,4 x 10 -2 ]. If 
we choose the diameter of the dimer subunit to be a = 5 nm, 
these concentrations correspond to 27-540 /iM. 

We note that Nguyen and coworkers l28tl recently devel- 
oped a model for capsids subunits in which subunit excluded 
volumes have a roughly trapezoidal shape. Interestingly, as- 
sembly kinetics predicted by their model are qualitatively sim- 
ilar to those of the model described above 12611 . except that 
they find insertion of the final subunit (to form a complete 
capsid) is rate limiting, and uphill in free energy for many 
sets of parameters. While insertion of the final subunit is also 
impeded by excluded volume for the present model, we find 
that once completed, capsids are stable and dissociation of a 
subunit is slow compared to assembly timescales. This result 
seems consistent with experimental observations that subunit 
exchange between completed P22 capsids and free subunits is 
characterized by longtimescales (days) compared to those for 
assembly (minutes) I691I7011 . 
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B. Modeling nanoparticle encapsidation 

We modify the empty capsid model by introducing a 
nanoparticle, or a rigid sphere, with radius Rc at a fixed posi- 
tion at the center of the simulation cell C. In addition to the 
pairwise interactions between subunits described above, sub- 
units interact with the nanoparticle via excluded volume in- 
teractions and attractive interactions. The potential energy of 
interaction uc(|R — C|) between a subunit at position R and 
the nanoparticle is a spherically symmetric shifted Lennard- 
Jones potential 

«c(r) = £cO(r c - f) [C(f/a) - C(r c /a)] (8) 

where r c = 2.5 is the cutoff distance. The strength of the at- 
tractive interaction is dictated by e c and f = r — (Rc — 0.5) 
is a shifted distance so that the attraction has its minimum 
value, — e c — u c (r c ) , when the center of the subunit and the 
surface of the nanoparticle are separated by (2 1 / 6 — 0.5)<r, 
and maintains the short range nature of the interactions con- 
sidered in Ref. [26]. This potential mimics core-subunit in- 
teractions that do not favor particular subunit orientations; we 
consider electrostatic interactions that depend on subunit ori- 
entations in a future work. The core-subunit interaction free 
energy, g c = e c — Ts a( j, includes an entropy penalty, s al j, for 
frozen degrees of freedom in the direction normal to the sur- 
face. We determine s a( j by calculating the partition function 
of an adsorbed subunit according to Eq. ©; the result varies 
weakly with the adsorption energy: 1 < s^/k^ < 2.4 for 
3<e c <12. 

Subunit reservoir. We represent a nanoparticle immersed in 
bulk solution without explicitly simulating thousands of sub- 
units by coupling dynamical simulations to a reservoir of sub- 
units at constant chemical potential. We divide the simulation 
box into a "main" region centered around the core, where or- 
dinary dynamics are performed, and an outer "bath" region 
where, in addition to ordinary dynamics, subunits are inserted 
or deleted in grand canonical Monte Carlo moves |p73tl . The 
main region is chosen to have a side length of L m = 15a, 
which is large enough that subunits cannot simultaneously in- 
teract with the nanoparticle and a subunit in the bath area, 
while the bath has a total width of L\, = 6a. Insertions and 
deletions are attempted with a frequency consistent with the 
diffusion limited rate for a spherical volume with a diameter 
of L m + L\,. As assembly proceeds, the concentration of free 
subunits is depleted and the bath chemical potential should be 
updated self-consistently. In this work, we consider assembly 
around a single nanoparticle in infinite dilution (i.e. Cc = 0) 
at system parameters for which little or no spontaneous as- 
sembly occurs away from cores, as we will see in the next 
section. Hence, the chemical potential remains constant. 

HI. RESULTS 

We have simulated assembly dynamics over ranges of sub- 
unit concentrations Cs, binding energies Eb, subunit speci- 
ficity parameters 6 m , and surface attraction strengths e c . All 
simulations use <fi m — tt and Rc = 1.2. 



A. The kinetics of core-controlled assembly 

In this section we present simple scaling arguments for the 
effect of core-subunit interactions on the kinetics of assembly 
and illustrate scaling in simulation assembly trajectories. We 
first concentrate on the average time to form a capsid, start- 
ing from unassembled subunits. As shown by Zlotnick and 
coworkers 11481 , 15011 . the assembly of empty capsids can often 
be broken into nucleation and elongation phases. We show 
elsewhere Jt^I that the average timescales of these phases for 
an individual capsid can be described by r — r nuc + r e i ong , 
with oc /4 W and » ci//(iV - n nuc ), where / 
is the subunit-subunit binding rate constant, c\ is the concen- 
tration of free subunits, n nuc is the number of subunits in the 
nucleus, and we assume c\ remains roughly constant during 
the assembly of an individual capsid. Because elongation re- 
quires N — n nuc assembly events, it introduces a minimum 
timescale for the overall assembly process, which is primar- 
ily responsible for the lag time in assembly kinetics reported 
in experiments J^, [l2l [5(3], theory l2ll 14811 . and simulations 
ll26l u§L l75tl . and results in a distribution of assembly times 
for an individual capsid that cannot be fit with a sum of pure 
exponential functions [74]. The observed assembly rate con- 
stant, /, can be considered an average quantity, since compu- 
tational models [26, 28] suggest that it varies for different in- 
termediates and decreases due to excluded volume constraints 
as assembly nears completion. In one model l28ll . insertion 
of the final subunit is slow compared to the rest of elongation 
and thus introduces a third timescale. Zandi, van der Schoot, 
and coworkers use continuum theory approaches to analyze 
nucleation ll76ll and capsid formation rates at long times OOll . 
Their finding that the total rate of capsid formation is pro- 
portional to cf at long times is consistent with the timescales 
given above for a single capsid if elongation dominates. 

Assembly rates on cores. As we will see from the simu- 
lations described below, the presence of cores modifies the 
nucleation and elongation timescales and introduces a new 
one, which describes the adsorption of subunits to the core 
surface. In this discussion we assume that cores are commen- 
surate with the size and geometry of capsids; we discuss the 
general case elsewhere. If there is no assembly of adsorbed 
subunits, the adsorption timescale is r a d = n s /fc a dCi, with k^ 
the adsorption rate constant, and the number of adsorbed sub- 
units at saturation, n s , can be calculated for Langmuir adsorp- 
tion 

n s = V c ci exp(-0g c )/ (1 + c\ exp(-j3g )) (9) 

where Vc ~ Na 3 is the volume available to adsorbed sub- 
units, (3 = 1/fcfiT is the inverse of the thermal energy, and 
g c is the surface-subunit free energy. Adsorption will usually 
be fast compared to assembly rates, which are slow compared 
to the diffusion limited rate for protein collisions 1 50] . Sim- 
ulation results demonstrate fast adsorption; Fig. [2] shows the 
number of adsorbed subunits, n a , as a function of time for 
several subunit concentrations. In each case, there is a rapid 
initial rise in the number of adsorbed subunits, which analysis 
of trajectories confirms is due to non-specific subunit adsorp- 
tion (i.e. without binding to other subunits). 
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FIG. 2: The time dependence of adsorbed subunits, n a , from simu- 
lations. Curves at increasing height correspond to reduced subunit 
concentrations of 10 3 Cs = 2.04, 4.07, 8.14, 20.4, 40.7, with a sur- 
face free energy of e c = 7 (<? c = —5.1), 8 m — 1, and Eb = 10. 
Simulation results are averaged over 60 independent trajectories. 

For fast adsorption, nucleation takes place at an effective 
surface concentration of c sur f = n s cr 3 /Vc with a timescale 

—core r — ritm firw 

T nuc - /surfC surf (10) 

where we define the surface assembly rate constant / sur f = 
fDc/D , with D and Dq the diffusion constants for free and 
adsorbed subunits, respectively, and we require that assembly 
rates are proportional to the frequency of subunit-subunit col- 
lisions [771] ■ In simulations for this work Dc = D because 
subunit friction is isotropic, but we will explore the effects of 
impeded surface diffusion elsewhere. Desorption of a nucle- 
ated intermediate is unlikely, since it would require breaking 
multiple subunit-subunit or subunit-core interactions. Assem- 
bly therefore leads to a positive flux of adsorbed subunits and 
elongation occurs at roughly the same concentration c sur f with 
a timescale 

Tetag ~(N- n s )/fc ad ci + (N — n nuc )// surf c surf (11) 

The first term on the right-hand side of equation ( [Til is the 
time for the remaining subunits to adsorb to the core surface 
while the second term accounts for the elongation reaction 
time. 

Eqs. [10] and QT| predict that cores enhance assembly rates 
by a factor 

r/r C ° re = T nuc + belong 

and that the relative timescales for nucleation and elongation 
can be manipulated by varying the surface-subunit free en- 
ergy, to yield regimes in which either nucleation or elonga- 
tion is rate limiting. Fig.[3]shows average assembly times for 
core encapsidation in simulations at varying surface energies 
4.1 < e c /k B T < 12; the upper and lower dashed lines iden- 
tify the scaling relations predicted by Eqs. ( fTTJb and ( fTTb for 
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FIG. 3: (a) Average assembly time, r, for simulations of core en- 
capsidation for varying surface energies plotted as a function of the 
surface concentration, csurf, described in the text. Simulated sur- 
face energies spanned e c £ [4.1, 12], giving surface free energies 
of — g c £ [2.7, 9.6]. The upper dashed line is a guide to the eye 
that indicates a nucleation dominated scaling of cf mf while the lower 
dashed line is a fit to the rightmost four data points with the form 
t = A + -B/c S urf to illustrate the elongation dominated scaling. Data 
points represent an average of 30 or more independent encapsidation 
trajectories, run at parameter values of Cs = 8 x 10 -3 , Eb = 10, 
and 9 m = 1. (b) Nine individual trajectories are shown for a surface 
energy in the nucleation regime, e c = 4.5 (c sur f = 0.15). 

nucleation and elongation dominated regimes, respectively, 
and we take n nuc = 5, although more data would be required 
to precisely estimate the nucleation size. In addition, the scal- 
ing relation is limited in range because the nucleation size can 
increase at very low e c . The nucleus usually corresponds to a 
small polygon in our simulations (see Fig. [5]), but the size and 
geometry of nuclei depend on system parameters. Time series 
of n a (the total number of adsorbed subunits) from individ- 
ual trajectories in the nucleation dominated regime are shown 
in Fig. [3J5 to illustrate the stochastic nature of the nucleation 
event. 

Cooperative assembly. For large surface free energies, 
Csurf ~ 1 , meaning that adsorption does not saturate until 



6 




0.5 1 1.5 

f/10 4 ? 




FIG. 4: The time dependences of the total number of adsorbed 
subunits, n a (solid lines), and the largest assembled cluster, n aS semb 
(dashed lines), reveal different assembly mechanisms. The subunit 
concentration is Cs = 8 x 1CP 3 , and the energy parameters are: (a) 
low surface energy and moderate binding energy, e c = 4.5, £b = 10; 
(b) high surface energy and low binding energy, e c = 12, £b = 7. 
The points labeled with (T) in (a) and (b) correspond to the structures 
shown in Figs. 5a and 5b, respectively. 



enough subunits have adsorbed to form a capsid, n s « N. 
Assembly in this regime can occur through a collective re- 
orientation of adsorbed subunits; at low subunit-subunit bind- 
ing energies (st>) this process typically takes place well af- 
ter adsorption has saturated and thus resembles the reverse 
micelle assembly mechanism suggested by McPherson l63ll . 
The time dependences of the number of adsorbed subunits, n a , 
and the number of assembled subunits, n assem b, are shown for 
representative trajectories that illustrate the different assembly 
pathways in Fig. [4] and structures from these trajectories are 
shown in Fig. [5] 



B. Packaging efficiencies 

In addition to enhancing and controlling rates of assembly, 
cores can increase assembly yields. As a measure of efficiency 
of assembly, we observe packaging efficiencies, which are de- 



(a) 




FIG. 5: Snapshots from the simulation trajectories shown in Fig.|4^ 
illustrate two assembly mechanisms: (a) sequential assembly at low 
surface energy, e c = 4.5, and (b) cooperative assembly at high sur- 
face enery, e c = 12. Snapshots from right to left correspond to in- 
creasing time and correspond to the triangles shown in (a) Fig. [4^ 
and (b) Fig. |4j;. The size of subunits is reduced to aid visibility, and 
subunit color indicates the number of complementary interactions: 
white, 0; green, 1; blue, 2; red, 3 or 4. All images of simulation 
structures in this work were generated with VMD [78]. 




6 8 10 12 14 



FIG. 6: The efficiency of encapsidation of a model nanoparticle 
is compared to the fraction of subunits in capsids due to sponta- 
neous assembly of empty capsids, shown with symbol (+). Packag- 
ing efficiencies, or the fraction of independent trajectories in which 
a nanoparticle was encapsidated by a perfect capsid, are shown for 
neutral subunits at e c — 7 (■) and e c = 12 (•), with a subunit 
concentration of Cs = 8 x 10~ 3 and a final observation time of 
tf = 48, 000. Spontaneous assembly results are shown at Cs = 0.11 
and tf = 600, 000. The specificity parameter 8 m — 0.5 for all points. 



fined as the fraction of independent trajectories for which a 
nanoparticle is encapsidated by a perfect capsid. A 'perfect 
capsid' is comprised of 30 subunits, each of which has the 
maximum number of bonds, 4. The variation of packaging 
efficiencies with binding energy is shown in Fig.|6]for an ob- 
servation time of tf = 48, 000, beyond which packaging effi- 
ciencies increase only slowly. 

Core control of assembly through heterogeneous nucle- 
ation. At the subunit concentration considered in Fig. [6] 
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Cs = 8 x 10~ 3 , spontaneous assembly into properly formed 
empty capsids is not observed for any of the binding energies 
considered, while cores are efficiently encapsidated over rela- 
tively wide range of e\>. This observation is consistent with ex- 
periments, which find that assembly occurs in the presence of 
nanoparticles below the critical subunit concentration (CSC) 
at which spontaneous assembly occurs [64]. Model nanopar- 
ticles enhance assembly because favorable core-subunit inter- 
actions lead to a high local concentration, c sur f, of 'adsorbed' 
subunits near core surfaces (see Eq. (0). Assembly and en- 
capsidation occur when the effective surface concentration ex- 
ceeds the CSC. Core encapsidation simulations were not car- 
ried out at binding energies of £b > 1 1 to ensure that there was 
no assembly in the bath. At higher binding energies, sponta- 
neous assembly is rapid and depletes the concentration of free 
subunits, and thus suppresses nanoparticle encapsidation and 
decreases packaging efficiencies. The competition between 
core-controlled and spontaneous assembly will be explored in 
a future work. 

Assembly on cores is robust in the sense that packaging ef- 
ficiencies remain near 100% over wide ranges of the subunit 
binding energy Eb and the surface attraction energy, e c . As 
a comparison, we consider independent simulations of empty 
capsid assembly (without model nanoparticles) at a subunit 
concentration of Cs = 0.11 for whichspontaneous assembly 
is relatively productive. As a measure of assembly effective- 
ness in empty capsid simulations, we define the 'packaging 
efficiency' as the fraction of subunits in complete capsids. At 
an observation time of tf = 6 x 10 5 , approximately 12 times 
longer than the observation time for core encapsidation, ef- 
ficient assembly occurs over the relatively narrow range of 
11.0 < £b < 13.5. Empty capsid assembly is thwarted by 
two forms of kinetic traps at higher values of subunit-subunit 
binding energies. If new assembly intermediates form too 
rapidly, the pool of free subunits becomes dep leted before 
most capsids finish assembling ifH [2lj |2r| |47t|48 t l49l |50ll , 
and malformed structures result when additional subunits bind 
more rapidly than strained bonds can anneal within a partial 
capsid i2q,l28il5lll52ll . Cores suppress the first of these traps 
by enabling rapid assembly well below the CSC, so that the 
number of nucleation sites is controlled by the concentration 
of cores even during rapid assembly. 

Core control of assembly through templating. The for- 
mation of mis-bonded configurations impedes assembly when 
subunit-subunit binding energies are large compared to k^T, 
because progression from these configurations to a properly 
formed capsid requires unbinding events and thus is charac- 
terized by large activation energies. An interior core could 
suppress this form of kinetic trap by acting as a template that 
directs assembly at all stages towards a morphology consistent 
with the low free energy capsid. We explore this capability by 
varying the subunit specificity parameter 9 m , which controls 
the intrinsic likelihood of subunits to form strained bonds (see 
section UTAl i. 

As shown in Fig. [7] efficient spontaneous assembly of 
empty capsids occurs over relatively narrow range of 0.5 < 
9 m < 1.0. At low values of 9 m , subunit binding rates are pro- 
hibitively slow because most collisions do not lead to bond 




FIG. 7: (a) Variation of assembly with specificity parameter, 9 m . 
Packaging efficiencies (•) are shown for Eb — 10, Cs = 8 x 1CP 3 
and t{ — 48, 000, while the fractions of subunits in empty capsids 
(+) are shown for e b = 12, Cs = 0.11, and. U = 240, 000. 



formation, while at larger values of 9 m strained bonds tend to 
be trapped within growing capsids. The presence of assem- 
blages with strained bonds and the lack of free subunits are 
illustrated by a snapshot from the end of a simulation with 
9 m = 2.0 in Fig. [8^; a snapshot at the end of a core simulation 
is shown in Fig. [8}? for comparison. 

Core controlled assembly, on the other hand, results in 
packaging efficiencies near 100% for specificity parameters 
as large as 9 m — 2.25 . Because the core is commensurate 
with the size of a perfect capsid, subunits are driven to bind 
with the correct local curvature at all stages of the assem- 
bly process. Note that the core simulations in Fig. [7] benefit 
from heterogeneous nucleation as well as templating; subunit 
binding rates and free energies increase with 9 m , and thus so 
do spontaneous nucleation rates. We determine that templat- 
ing becomes increasingly important as 9 m rises because empty 
capsid simulations with 9 m > 1.5 yield a significant fraction 
of malformed structures. Simulations to be presented at a fu- 
ture work, in which there is heterogeneous nucleation but no 
templating, show a similar degree of sensitivity to 9 m . 



IV. DISCUSSION 

In this work we present scaling arguments and simula- 
tions that describe the assembly of capsid protein subunits 
around rigid cores. The kinetics and efficiency of assembly 
are predicted as functions of subunit-subunit binding energies, 
subunit concentrations and subunit-core interaction strengths. 
We find that assembly pathways depend sensitively on the 
strength of subunit-core interactions. For weak attractions, 
capsid formation requires assembly of adsorbed subunits into 
a stable intermediate, followed by sequential adsorption and 
assembly of individual subunits. Strong interactions induce 
rapid adsorption resulting in nearly complete, but disordered, 
monolayers of subunits, followed by cooperative subunit re- 
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FIG. 8: (a) Snapshot of part of the simulation box at the end of an 
empty capsid simulation with 9 m — 2 for the parameters in Fig. [7] 
Some capsids are closed but strained, as indicated by the presence of 
hexagons or squares, (b) Snapshot at the end of the core encapsida- 
tion trajectory shown in Fig.|5^. 



ordering to form capsids. Variations in assembly mechanisms 
are revealed by the time dependence of adsorbed subunits, as 
shown in Fig. [2] and by the scaling of assembly times with 
subunit concentration, as shown in Fig. [3] 

Model predictions for the variation of assembly timescales 
and packaging efficiencies with subunit concentration, Cs, 
and core-subunit interaction strengths, g c , can be verified in 
virus-like particle assembly experiments. The parameter g c 
can be related to the experimentally controlled functionalized 
surface charge density. Protein adsorption and assembly ki- 
netics can be monitored with time resolved light scattering 
[64], and packaging efficiencies can be determined from TEM 
micrographs 13311 . Successful validation of model predictions 
would provide strong evidence for the control of assembly 
mechanisms by a templating component, which prior works 
suggest is important for understanding viral assembly in vivo 

Implications for desi gnin g and understanding asse mbly re- 
actions. Experiments 11121 pill and models 112 ll |26j |28j |48j 
l49l I75I1 show that subunits can spontaneously assemble into 
low free energy ordered states with high yield and selectiv- 
ity, but that effective assembly is limited to optimal ranges 
of the forces that drive assembly. For non-optimal interac- 
tions, assembly is either not thermodynamically favored, or 



thwarted by long-lived disordered states. These limitations 
of spontaneous assembly have shaped the evolution of assem- 
bling components in biological systems, and similarly con- 
strain the development of assembled nanostructured materials. 
Our results suggest that cooperative interactions between dis- 
parate assembling components offer the potential to circum- 
vent some limitations of spontaneous assembly, particularly 
through heterogeneous nucleation and templating. 

The novel mechanisms and capabilities of multicomponent 
assembly introduce new considerations for the design of as- 
sembly processes. For example, Zlotnick and coworkers ll50ll 
show that a slow nucleation step in the spontaneous assembly 
of empty capsids can suppress kinetic traps. This condition 
is met for some parameter values our simulation model - ini- 
tial assembly steps can be slow in comparison to later ones 
because subunits in small intermediates have few bonds. As 
subunit-subunit binding energies are decreased these initial 
steps become even slower; however, the capsid products are 
less thermodynamically stable and subsequent "elongation" 
rates can also decrease, resulting in increased assembly times. 
Introducing cores to the simulations provides an independent 
means to control the formation of assembly nuclei while still 
enabling fast elongation kinetics, and thereby promotes rapid 
assembly with high yield and selectivity. 

Controling assembly kinetics by changing the properties of 
a template component will be useful for designing synthetic 
or biomimetic assembly reactions for which it is impractical 
to change the molecular structure of subunits or environmen- 
tal conditions. For example, protein-protein interactions, and 
hence the critical subunit concentration (CSC), can be con- 
trolled in virus-like particle experiments by varying the salt 
concentration or pH lfl2l l29i RtIi . Capsid proteins denature, 
however, if these parameters are changed too far from physi- 
ological conditions. Our results suggest that varying the func- 
tionalized surface charge density on nanoparticles enables in- 
dependent control over the CSC and packaging efficiencies. 
Similarly, viruses have limited capability to control the cel- 
lular environments in which they replicate, and amino acids 
at capsid protein-protein interfaces are highly conserved, per- 
haps in part because the need for capsid dissociation upon in- 
fecting a new cell constrains these interactions. Interactions 
between capsid proteins and the viral genome or host cellu- 
lar compenents may provide important alternative avenues to 
promote and control assembly, and hence viral replication. 



V. OUTLOOK 

Core-assembly geometry incompatibility. Our models de- 
scribe encapsidation of cores with shapes and sizes commen- 
surate with the low free energy capsid product. The simulation 
results demonstrate that the influence of core curvature on lo- 
cal subunit-subunit bonding configurations can dramatically 
influence global capsid morphologies (see Fig. |7). Core cur- 
vature that is inconsistent with the lowest free energy subunit- 
subunit bonding configurations will introduce frustration and 
thus may limit the robustness of assembly. Experimental ob- 
servations that solid cores l33ll and nucleic acid cores ll79ll 
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with different sizes promote assembly of different capsid mor- 
phologies demonstrate that frustration is an important consid- 
eration in biological and nanostructured assembly processes. 
Comparison of our current model with experimental results 
will identify frustration and future work will explicitly ad- 
dress geometrical frustration. Additionally, physio-adsorption 
of subunits on core surfaces could impede lateral diffusion and 
thereby promote kinetic traps. However, simulations in which 
adsorbed subunits had friction constants increased by a fac- 
tor of 100 in directions tangential to the core surface demon- 
strated only a small increase in propensity for kinetic traps, 
although net assembly rates were slower. 

Fluctuating cores. Cores comprised of nucleic acids 
or other macromolecules with dynamic configurations can 
change size and shape during encapsidation. Although the 
timescale arguments presented in Eqs.[lO]and[TT]can be gen- 
eralized to include an additional timescale that represents core 
dynamics, the possibility of additional forms of kinetic traps 
due to fluctuating core configurations should be explored. 

VI. CONCLUSIONS 

In summary, we have presented simulations and scaling ar- 
guments that describe the assembly of solubilized subunits 



around rigid cores. These models mimic the dynamical as- 
sembly of viral capsid proteins around functionalized inor- 
ganic nanoparticles, but are general enough for broad appli- 
cability in describing the assembly of biological or nanos- 
tructured materials around templates. We find that template 
properties can dramatically influence assembly timescales and 
mechanisms, as evidenced by the prediction of a novel assem- 
bly mechanism not seen during the assembly of empty cap- 
sids. These conclusions may be significant for understanding 
the role of nucleic acids in the viral assembly and for design- 
ing nanomaterials or drug delivery vehicles to interact with 
cargo molecules. 
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